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ABSTRACT 

We present a systematic study of the evolution of intermediate- and low-mass 
^ • X-ray binaries consisting of an accreting neutron star of mass 1.0 — I.SMq and 

a donor star of mass 1.0 — G.OMq. In our calculations we take into account 
physical processes such as unstable disk accretion, radio ejection, bump-induced 
detachment, and outflow from the L2 point. Comparing the calculated results 
with the observations of binary radio pulsars, we report the following results. 
(1) The allowed parameter space for forming binary pulsars in the initial orbital 
period - donor mass plane increases with increasing neutron star mass. This may 
help explain why some MSPs with orbital periods longer than ~ 60 days seem to 
have less massive white dwarfs than expected. Alternatively, some of these wide 
binary pulsars may be formed through mass transfer driven by planet/brown 
dwarf- involved common envelope evolution. (2) Some of the pulsars in compact 
binaries might have evolved from intermediate-mass X-ray binaries with anoma- 
lous magnetic braking. (3) The equilibrium spin periods of neutron stars in 
low-mass X-ray binaries are in general shorter than the observed spin periods 
of binary pulsars by more than one order of magnitude, suggesting that either 
the simple equilibrium spin model does not apply, or there are other mecha- 
nisms/processes spinning down the neutron stars. 

Subject headings: stars: millisecond pulsars — stars: evolution — MSPs: general 



1. Introduction 



Millisecond pulsars (MSPs) are neutron stars (NSs) characterized by short spin periods 
(-Pspin < 20 ms) and weak surface magnetic fields {B < 10^ G), which are often found in 
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binaries with a white dwarf (WD) companion. It is generally agreed that MSPs are old 
NSs, recycled by accretion of mass and angular momentum from the donor stars through 
Roche-lobe overflow (RLOF) during the previous low-mass X-ray binary (LMXB) evolution. 
The NS was spun up, due to mass accretion, into a MSP, while the donor evolved to be a He 



or CO WD (see e.g., iBhattacharya fc van den Heuvellll99ll : iTauris fc van den Heuvelll2006l . 
for reviews). 

The stability of mass transfer in an X-ray binary depends on the ratio of the masses 
of the donor star and the NS, and the initial orbital period of the system. Tradition- 
ally, it was thought that, if the mass ratio is large enough (> 1.5) , the mass transfer is 



likely to be uns t able, resulting in a common envelope (CE) evolution ( iPaczynski et al.lll976 



WebbinkI Il984j : llben fc Livid Il993l ). where the NS spirals into the env elope of the donor 



on a very short timescale (< 10^ yr). More recent investigations (e.g. 



Podsiadlowski fc RappaportI l2000l : iKolb et al.l l2000l : iPodsiadlowski et al. 



auris et al 



20021 : 



200 



Pfahl et al. 



20031 ) show that X-ray binary systems with intermediate-mass (up to ~ 5Mq) donor stars 
(i.e., IMXBs) can avoid the spiral-in phase and experience rapid mass transfer on a thermal 
timescale, successfully evolve to become LMXBs. 

For LMXBs, a CE evolution is even unavoidable if the separation of binary components 
is large enough, so that the donor reaches the asymptotic giant branch (AGB) phase with a 
deep convective envelope before RLOF (i.e. Case C RLOF). Stable mass transfer in LMXBs 
usually occurs on a timescale ~ 10® — 10^" yr when the donor star is on the main sequence or 
(sub)giant branch at the onset of RLOF (i.e.. C ase A or B RLOF), a r id for r ning M SPs seems 
to be feasible in this process. It was found by iPylyser fc Savonijd (Il988l . Il989l ) that there 
exists a critical bifurcation orbital period (Pbif), the initial orbital period that separates the 
formation of converging LMXBs (which evolve with decreasing orbital periods until the donor 
star becomes degenerate or an ultra-compact binary is formed) from diverging LMXBs. The 
value of Pbif is found to be ~ 1 d, but depends heavily on the processes of tidal interactions 



and the mechanisms and efficiency of orbital angular momentum loss (lErgma et al.l Il998 



Podsiadlowski et al.l 120021 : Ivan der Sluvs. Verbunt. fc Pold 120051 : iMa fc Li 1120091 ). 



The final products of LMXB evolution are binary pulsars. The distributions of the 
spin periods of the pulsars, the or bital pe r iods a nd the WD masses can be used to testify 
the models of I/LMXB evolution. iDeloyd ( 120081 ) compared the theoretical expectations of 
I/LMXB evolution to the populations of Galactic binary pulsars. He showed that a significant 
population of binary pulsars with 1 d < Porb ^ 100 d are generally consistent with being the 
descendants of long-period LMXBs or IMXBs. However, there remain quite a few unresolved 
puzzles. For example, binary pulsars with Porb ^ 60 d seem to have WD compani o ns les s 
massive than predicted by theory, as pointed out previously by iTauris &: Savonijd (119991 ). 



- 3 - 



and those with 0.1 d < Porb ^ 1 d are inconsistent with any I/LMXB evolution. 



In the previous studies on I/LMXB evolution, a canonical NS (of mass ~ 1.3 — 1.4 Mq 
was usually adopted. This seems to be supported by the finding that the NS mass distri- 
bution is consistent with a narrow Gaussian at 1.35 ± 0.04 (IThorsett fc Chakrabartv 



1999 ). However, both observations fsee IZhang et al.l 



2010l : ISchwab. Podsiadlows. 



and theories (INomoto 



d. fc Rappapor 



1984 



Timmes et al 



201C; 



1996 



2011 



Ozel et al. 



Heger et al 



Kiziltan. Kottas. &: Thorsett 



20121, and references therein) 



20031 : IWooslev fc Jankall2005h 



suggest that the initial masses of NSs may occupy a large range, probably originating from 
two different mechanisms of forming NS: iron-core collapse supernovae and electron-capture 
supernovae. NSs in high-mass X-ray binaries (HMXBs) have experienced very little accretion 
because of their young ages, so their masses should be very close to those at birth. The mea- 



sured masses of NSs in HMXBs range from l.OGl^ ln M^ for SMC X-1 fivan der Meeret al.l 



2007h to 1.8 6 ±0.1 6 M© for Vela X-1 jBarziv et allboOll : lOuaintrell e7aDl2003[ l. Recentlv 
Rawls et al.l (120111 ) present an improved method for determining the mass of NSs in eclips- 
ing X-ray pulsar binaries and apply it to six systems. They find that the NS masses range 
from 0.87 ± 0.07 Mq (eccentric orbit) or 1.00 ± 0.10 (circular orbit) for 4U1538-52 to 
1.77 ± 0.08 M0 for Vela X— 1. So in a proper investigati on, the infiuence of th e NS m asses 
should be included. Recent evolutionary calculations by iDe Vito fc Benvenutd (120 lOl ) have 
provided evidence that the evolution of I/LMXBs depends upon the NS mass. In this work 
we perform systematic calculations of I/LMXB evolution and discuss the properties of the 
produced binary and millisecond pulsars (BMSPs), taking into account different initial NS 
masses. We adopt the initial donor masses to be 1.0 — 6. OM©, and two different initial masses 
(1.0 and 1.8 Mq) for the NS. 

During the mass transfer processes, part of the transferred mass from the donor star 
may escape from the binary system, carrying away the orbital angular momentum. The 
formation of BMSPs is closely related to the mechanisms of mass and angular momentum 
loss, and the stability of ra ass transfer is also dependent on the angular mementum loss 
rate (ISoberman et al.l 119971 ) . Generally, it is assumed that all the mass transferred from 
the donor will accrete onto the NS unless in a super-Eddington mass transfer phase, during 
which the NS will accrete at the Eddington rate (MEdd ~ 1-5 x lO'^Mgyr"-^ for a 1.4 Mq 
NS), and the residual mass escapes from the binary system, carrying the NS's specific orbital 
angular momentum. This is the so-called "isotropic reemission model" . Actually, there may 
also be mass loss even during the sub-Eddington mass transfer phase. For example, the 
accretion disk can be come thermally and viscously unstable when the orbital period is larger 
than a critical value ( van Paradijs 1996 : Dubus et al.lll999l ). leading to limit cycle behavior 
of the mass transfer rate. The NS can accrete mass only during outbursts, while most 
matter may be ejected out of the binary systems during quiescence by the radiation and 
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magnetic pressure of the rapidly rotating NS (iRuderman et al.l Il989l : iBurderi et al.ll2002l ). 
In our calculation s, mass loss due to s uper-Eddington mass transfer, r adio ejection caused b y 
an unstable disk (IBurderi et al.ll2002l ) or bump-induced detachment (jP'Antona et al.ll2006l ). 
and in some cases outflow from the L2 Lagrangian point are included. 

This paper is organized as follows. In Section 2 we describe the stellar evolution code 
and the binary model used in this paper. We present the calculated results in Section 3 and 
discuss their possible applications in the formation of BMSPs in Section 4, and summarize 
in Section 5. 



2. Binary Evolutionary Calculations 

2.1. The Stellar Evolution Code 

All calcul ations we r e carr i ed ou t with an updated version o f the stellar evolu tion code 
developed by lEedetonl fll97ll . Il972h (see also Iflan et al.l Il994j : IPols et al.lll995h . We set 
initial solar chemical compositions (i.e., X = 0.7, Y = 0.28, and Z = 0.02) for the donor 
star, and take the ratio of the mixing length to the pressure scale height to be 2.0, and the 
general convective overshooting parameter to be 0.12. In our calculations, we have considered 
a number of binary interactions in order to follow the details of the mass transfer process, 
including orbital angular momentum loss due to gravitational radiation (GR), magnetic 
braking (MB), mass loss, and the effect of disk instability. 



2.2. The Input Physics 



Each of the binaries initially consists of an NS of mass Mi (= 1.0 or 1.8 Mq) and a 
zero-age main-sequence (ZAMS) donor of mass M2 (~ 1.0 — 6.0M (7)). The efective radius of 
the Roche lobe Rl oi the donor is given by the following formula (lEggletoru 119831 ) . 



Ri/a 



0.49g 



-2/3 



0.6g-2/3 -Mn(l + g-V3) ' 



where a is the orbital separation of the binary, and q = M2/M1 is the mass rat io of the bin ary 
components. As usual we assume that tides keep the binary orbit circular ( lKinglll988l ). so 
that the orbital angular momentum is 



(2) 
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where M = Mi + M2, and = \fGM~[c? is the orbital angular velocity. Logarithmic 
differentiation of Eq. (2) with time gives the rate of change in the orbital separation 



Ml 



M2 M 
2 =. -I 

M2 M' 



a Jorb Ml 

where the total changing rate in the orbital angular momentum is 



(3) 
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The three terms on the right-hand-side of Eq. (4) represent angular momentum loss due to 
GR, MB, and mass loss, respectively. The r ate of angular momentu r n loss due to G R is 
calculated according to the standard formula fiLandau fc Lifshitall959t lFaulknerlll97l[ ) 



J. 



GR 



J. 



orb 



M1M2M 



(5) 



where G and c are the gravitational c onstant and the speed of light, respectively. The 
prescription of IVerbunt fc ZwaanI (Il98ll ) is adopted to calculate the angular momentum loss 
due to MB, 

Jmb 



J, 



orb 



-3.8 X 10-3°^^ s-\ 



(6) 



where Ro is the radius of the donor. 



For IMXBs and wide LMXBs, the mass transfer rate IM2I (note that M2 < 0) can be 
larger than the Eddington accretion rate Mscid- Thus the mass loss rate from the binary 
system is M = M2 + ^Edd- Here, we adopt the isotropic reemission model, assuming that 
the extra material leaves the binary in the form of isotropic wind from the NS, carrying 
off the NS's specific orbital angular momentum ji = (M2/M1M) Jorb- The related angular 
momentum loss rate can be derived to be 



Jr 



ml,l 



Mo 



Jorb MiM 



(M2 + MEdd) 



(7) 



In some cases part of the material lost from the donor star may escape the system 
through the L2 Lagrangian point rather accrete onto the NS. Assuming that a fraction 6 of the 



matter flow escapes from the binary with the specific orbital angular mo mentum 72 



the angular momentum loss rate due to the L2 point outflow is given by (IShao &: Lill2012l ) 



J 



ml,2 



J 



orb 



gj, MM2 
a2 MiM2^ 
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where is the distance between the mass center of binary and the L2 point. 



The transferred material from the donor will form an accretion disk surrounding the 
NS. If the effective temperature in the accretion disk is below ~ 6500 K (the hydrogen 
i onization tem perature), the accretion disk is likely to be thermally and viscously unstable 
( Lasotall2001l) . In LMXBs irradiation from the NS may help stab ilize the disk to some extent 

J 



van Paradijs 1996 : King et al. 1997 : Lasota 2001 : Ritterl 2008 ). The critical mass transfer 



rate for the disk instability is given by (IDubus et al.lll999l ) 



3.2 X 10 
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(9) 



If the mass transfer rate is less than M^, the X-ray binary is assumed to become transient, 
experiencing short outbursts separated by long quiescent intervals. Material accumulates in 
the disk during the quiescent pha se while the NS a ccretes mostly during outbursts. Define 
the duty cycle d (~ 0.001 — 0.1, see King et al.ll2003l ) to be the ratio of the outburst timescale 
to the recurrence time, the accretion rate during outbursts can be estimated as ~ —M2/d 
. Of course, the mass accretion rate in this case is also limited by the Eddington accretion 
rate. 



3. Results of I/LMXB Evolution Calculations 
3.1. The Initial M2 — Porb Parameter Space for Successful Evolution 



In Fig. 1 we outline the results of our calculations b y showing the fate o f the I/LMXB 
evolution in the initial M2 — Porb diagram, similar as in iTauris et al.l (120001 ). The left and 
right panels correspond to 1.0 Mq and 1.8 Mq NS, respectively. It is seen that the allowed 
space for successful evolution into binary pulsars (i.e., without CE evolution) is larger for 
Ml = 1.8 Mq than for Mi = 1.0 Mq, since the lower mass ratio in the former case can stabilize 
mass transfer during the IMXB evolution. When Mi = 1.0 Mq, systems with a low mass 
companion (M2 < I.SMq) can form binary radio pulsars with a He/CO WD companion if 
the initial orbital period is < 1000 d, while the lower limit of M2 increases to ~ 2.2 Mq when 
Ml = 1.8 Mq. IMXBs can avoid spiral-in and CE evolution when 1.5Mq < M2 < 3.3Mq and 
2 d < Porb < 10 d for Ml = 1.0 Mq, and 2.4Mq < M2 < 5.5Mq and 2 d < Poa < 40 d for 
Ml = 1.8 Mq, respectively. In both cases, if the initial Porb is too short, the binary systems 
will experience eith er a CE phase, or become X-ray binaries with degenerate hydrogen stars 
(lErgma et al.l Il998l ). possibly forming ultra-compact binaries. On the other hand, if the 
initial Porb is too long, the donor will develop a deep convective envelope at the onset of 



RLOF, and a runaway mass transfer is initiated, leading to CE evolution. In iTauris et al. 
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(120001) systems with Mi = 1.3Mq, 2Mq < M2 < SM© and 1 d < Porb < 20 d can survive 
unstable mass transfer. These systems are right between the hmiting cases presented in 
Fig. 1. 

To see in more detail how the NS mass affects the binary evolution, we show in Fig. 2 
the evolutionary paths of two IMXBs with the same donor mass M2 = AMq and orbital 
period Porb = 10 d, but different NS mass (1.0 Mq and 1.8 Mq). At the age of ~ 175.1 Myr, 
the donor star starts to transfer mass via RLOF. In the upper two panels with a 1.0 Mq NS, 
mass transfer occurs rapidly, rising up to > lO'^M© yr~^ within ~ 10^ yr. Meanwhile, the 
orbital period reduces to < 4 d. This dynamically unstable mass transfer will result in CE 
evolution and probably merging of the NS and the He core of the donor. In the lower two 
panels with a 1.8 Mq NS, mass transfer initially proceeds at a rate of ~ 3 x 1O~^M0 yr~^ for 
~ 10^ yr, during which most of the hydrogen-rich envelope (~ ^ Mq) is removed from the 
donor star. The mass ratio inverts during this phase, and subsequently the mass transfer 
rate decreases to be several lO'^M© yr~^ for ~ 10^ yr. During this time the NS accretes 
about O.OI2M0 mass, and this is the primary process to spin up the NS. The final product 
is a binary consisting of a 1.815 Mq (mildly) recycled NS and a ~ 0.6Mq CO WD in an 
orbital period of ~ 36 d. 



3.2. The P/X - MwD Diagram 

Figure 3 shows the calculated correlation between the final orbital period P^b"^' 
the mass of the WD (the remnant of the donor) M-wd, for two different initial NS masses, 
I.OM0 (left panel) and I.8M0 (right panel). As mentioned above, low-mass donor stars 
(M2 < I.3M0 or < 2.2M0 for I.OM0 or I.8M0 NS, respectively) may evolve to be He 
WDs. However, when the mass of He is accumulated to exceed ~ 0.4 — O.5M0 in the core 
of the donor. He flash will occur, giving rise to the formation of a CO core. The orbital 
period will reach > 10^ days. Intermediate-mass donors can avoid He flash due to their 
higher temperature, forming CO WDs with the mass Mwd ^ 0.33 Mq. Their P/X"' ~ ^wd 
distribution deviates from that of the low-mass branch obviously. 



3.3. Radio Ejection during the "Bump-related" Detachment? 



When its hydrogen shell reaches a discontinuity in the hydrogen content at the time of 
the first dredge-up, a low-mass star will suffer a ten aporary contract i on, th us producing a 
"bump" in the luminosity function of the red giants. iD'Antona et al.l (120061 ) suggested that 
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the orbital period gap (~ 20 — 60 days) of BMSPs may be related to the bump-induced 
detachment of the donor star from its RL. If the NS has already accreted sufficient mass, it 
may turn on as a MSP during the detachment phase. Material transferred from the donor, 
once it expands again to refill its Roche lobe, may b e inhibited by the pulsar's radiation 



press ure and ejected at the inner lagrangian point Li ( jRuderman et al.lll989l : iBurderi et al. 



2OO2I ). so that no further mass accretion would occur. Similarly, this s o-called "radio eject ion" 
process may also occur in LMXBs with an unstable accretion disk ( iBurderi et al.ll2002l ). In 
both cases the mass loss will infiuence the orbital evolution of the binary and the spin 
evolution of the NS. iD'Antona et al.l (|2006[ ) considered the radio ejection only during the 
bump-related detachment. Our calculations show that the disk instability generally appears 
earlier than the bump phase. Two examples are presented in Fig. 4. In the upper two panels, 
the binary system initially consists of a 1.0 Mq NS and a 1.0 Mq ZAMS donor star in an 
orbital period of 5.0 d. RLOF initiates at the age ~ 1.186 x 10^" yr, and the bump occurs 
~ 7x lO'' yr later, at which the donor mass has decreased to be ~ O.52M0. However, the disk 
becomes thermally unstable when the donor mass is ~ O.OM©, much earlier than the bump 
phase. If the NS's spin was accelerated during the prior mass transfer, the accumulated mass 
in the disk during quiescence will be ejected out of the binary by the radiation pressure of 
the rapidly spinning NS. We find a similar situation in the case of a binary containing a 
1.8 Mq NS and a 1.5 Mq donor star shown in the lower two panels. Compared with the 
effect of unstable disk accretion, it seems that the bump-related detachment might play a 
less important role in the evolution of LMXBs. 

We finally note that the process of pulsar-driven mass ejection is highly uncert ain, since 
its condition and efficiency depend on several unknown parameters (IFu fc Lill201ll ). such as 
the value of the equilibrium period that a NS will reach during the mass transfer. However, 
as seen below, there is controversy on the estimate of the equilibrium period for NSs in 
LMXBs, and one should be cautious when considering the effect of radio ejection on the 
binary evolution. 



3.4. The effect of outflow from the L2 point 

In the binary evolution the stability of mass transfer strongly depends on mass loss 
and related angular momentum loss. Although we have taken into account various ways 
for mass loss, including super-Eddington mass transfer, unstable mass transfer due to disk 
instability, and mass ejection due to the propeller effect and the pulsar's radiation, the 
detail processes are complicated and uncertain, and simplified treatment might not reveal 
the realistic situation. In the following, we show an example of how outfiows from the 
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binary systems through the L2 point influence the evolution of mass transfer, and change 
the conditions of forming BMSPs. 



MSP J1614-2230 is a 3.15 ms pulsar of mass 1.97 ± 0.04 M© f lDemorest 1120101 ). The 
mass of its WD com panion (0 . 500 =b 0.006 M(7)) and th e orbit al period (8.7 d) were also 
measured accurately. iLin et al.l ( 1201 ll ) and iTauris et al.l (120 111 ) systematically investigated 
the formation channels of this pulsar from an IMXB. They showed that NS this massive are 
not easy to produce in spite of the initially high r nass of the donor st ar, unless they were 
already born as a relatively massive NS. However, iTauris et al.l ( 1201 ll ) found that, for the 
system with a 1.8 Mq NS, the final orbital period will be always larger than 10 d, inconsistent 
with the observation of the pulsar. This conclusion will not hold i f we let a small fra ction 
6 of the transferred mass leave the system from the L2 point (e.g., iBailyn et al.lll989l ). As 



an illustration, in Fig. 5 we compare the evolution of an IMXB containing a 1.8 Mq NS and 
a 4.5 Mq companion star in an orbit of 2.4 d, with 6 = (upper two panels) and 6 = 0.04 
(lower two panels), respectively. In the former case the primary angular momentum loss 
mechanism is the isotropic reemission around the NS during the super-Eddington accretion 
phase. The final system consists of a 2 Mq NS and 0.5 Mq CO WD, with a orbital period 
15 d. In the latter case angular momentum loss due to the L2 point outflow plays a role 
when the system evolves into an LMXB, and reduces the orbital period to ~ 8.7 days. 



4. Comparison with BMSPs 
4.1. The Porb — relation for low-mass BMSPs 



In wide LMXBs, the donor will climb to the red-giant branch (RGB) in the HR diagram 
before RLOF. For low-mass stars (< 2.3M0) on the RGB, there is a well known relationship 
between the mass of the degenerate He core and the radius of the g iant star, which is 



almo s t entirely independen t of the mass of the hydrogen-rich envelope (|Refsdal fc Weigert 



I971I : IWebbink et al.lll983l ). Based on this relationship, a specific correlation between the 



orbital period Pprh and the W D mass Mwd is obtained (jjoss et al.l Il987l : iRappaport et al. 



I995I : ITauris &: Savonijelll999l ). Comparison with the observations shows that a significant 
population of BMSPs with He WD companion is generally consistent with this Porb — 
relation. However, there seems to be a systematic deviation from the correlati on for pulsar s 



)r p 

with Porb ^ 60 d, which seem to have WD companions lighter than expected (lTaurislll996l ). 



Both systematic small values of the orbital inclination i and large NS mass can increase 
Mwd for the given observed mass functions. Since there does not seem to be any obser- 
vational selection effect favoring small inclination angle i (iTaurisI Il996l ). we first examine 
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whether the Porb — Mwd correlatio n can be acc o unted for if the long-p eriod BMSPs have 



CO W D companions. For example, IStairs et al.l ( 120051 ) noticed that the iTauris fc Savonije 



( I999I ) Porb — My/D relation is incompatible at the 99.5% level with a uniform distribution 
of cosi if the pulsar masses are drawn from a Gaussian distribution centered on 1.35M0 
with width O.O4M0, and better agreement with uniformity in cosi can be reached if the 
pulsar masses are large on average (e.g., 1.75 ± O.O4M0). An extreme example is PSR 
B0820+02, which has a .6 Mq CO WD companion in a very wide orbit with Porb — 1232 d 
jKoester fc ReiinerJgOOgi) . However, in such wide binaries it is difficult for the NS to accrete 



enough matter (ILi &: Wanglll998l : iTauris &: Savonijelll999l ). and the NS must be born heavy 



if this interpretation is correct. In Fig. 6 we compare the relations between Porb and M-wd 
in the cases that the NS has an initial mass of 1.0 Mq (left panel) and 1.8 Mq (right panel). 
Also plotted are binary pulsars with measured Porb and Mwd (90% probability mass range 
for randomly oriented orbits) for a fixed NS mass of 1.2 and 2.0 Mq, respectively (i.e., in 
each case the NS is assumed to have accreted O.2M0 during the mass transfer). It is seen 
that some binary pulsars with Porb > 100 d can fairly match the relation if they have massive 
NSs (~ 2M0) and heavy CO WDs (although in some cases small orbital inclination angles 
may be required), while for those with Porb < 20 d, statistically lighter NSs (~ I.2M0) seem 
to follow the relation better. 

It is expected that the NS has to accrete at least a few ~ 0.01 Mq mass t o reach a mil- 



liseco nd period, but this is difficult to achieve for NSs in wide binaries (see also iLiu &: Chen 



20111 ). because the mass transfer rate (which increases with increasing orbital period) is 
likely to be super-Eddington, and the accretion disk is likely to be unstable. In Fig. 7 we 
show the mass transfer rate IM2I, and the accreted mass AMi of the NS as a function of the 
final orbital period Porb- Note that here we plot the mass transfer rate only for the stable 
mass transfer phase - when the accretion disk becomes unstable, we take it to be the critical 
value Mcr at the onset of disk instability. Hence both the mass accretion rate and AMi are 
limited by the value of Mcr'- when IM2I is less than the M^r, unstable disk accretion occurs, 
and the NS is assumed to accrete mass only during outbursts, and part of the matter will 
be ejected out of the binary systems if \M2/d\ is super-Eddington. 

The left panel of Fig. 7 shows the results for a 1.0 Mq NS with a companion star of initial 
mass M2 = 1.0, 2.0 and 3.0 Mq, respectively. Generally more massive donor stars result in 
higher mass transfer rate, which also increases with longer orbital period Porb- It is noted that 
(1) IM2I is always > O.lMEdd, and (2) AMi < O.3M0, which decreases with Porb, because in 
wider systems the mass accretion rates (with both stable and unstable accretion disks) are 
more likely to be super-Eddington. The right panel is for the systems with a I.8M0 NS. In 
the case of light companion star (~ I.OM0), the initial orbital period should be larger than 
~ 1 d, so that the critical mass transfer rate Mcr is always larger than the mass transfer 
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rate, since Met increases during the mass transfer pro cess (with increasing orbit period), 
while the mass transfer rate IM2I decreases all the way ( iWebbink et al.lll983l ). The accreted 
mass AMi by the NS is thus significantly lower than those with more massive donors. So 
only results with 2.0 and 3.0 Mq donor stars are presented, in which AMi < O.SSMq, and 
M2 > O.SMEdd. 

Figure 7 shows that in general M2 ^ O.lMEdd) and AMi < 0.3Mq. The distribution of 
the mass transfer rate is conside rably higher than that (M2 ~ 10~^^ — 1O~^M0 yr~^) obtained 
by iPodsiadlowski et al.l (120021 ) . but more compatible with the observations of persistent 
LMXBs. The main reason is that we have tak en into account the effect of disk instability. 
The AMi distribution is roughly in line with Liu fc Chen (j2011 ). who found that AMi 
is generally less than 0.6 Mq in their calculations for systems with a 1.4 Mq NS and a 
1.0 — 2.OM0 donor star. Moreover, we find that systems with initially massive NSs (1.8 Mq) 
may accrete enough mass to evolve into BMSPs with 60 d < Porb ^ 200 d, , while light NSs in 
wide binaries are more likely to be partially recycled. This seems to support our conjecture 
that some of the wide BMSPs might be born massive. Hopefully accurate measurements of 
both the pulsar and the WD masses will help settle this issue. 



Lohmer et al.l (120051 ) performed timing observations of the BMSP J1640+2224 (with an 



orbital period of 175 d), and constrained the WD mass to be 0.151q;o5 Mq (Icr uncertainties), 
which indicates that the companion is very likely to be a low-mass He WD. In this case 
the massive NS + CO WD model obviously does not work, and one has to explore other 
possibilities. Evaporation of the companion star from a wind of relativistic particles after 
the pulsar turns on may decrease the companion mass significantly, but it is unlik ely for 
wide BMSPs, since the evaporation timescale would be longer than the Hubble time (ITauris 
1996h . 



Here we suggest another possible solution. It is interesting to note that solar-type stars 
are usual ly found to be surro unded by sub-stellar companions (usually planets and/or brown 
dwarfs) (jCassan et al.ll2012l ). One may expect that in some relatively wide LMXBs the 
companion star had possessed substellar companion(s) in close orbits like "hot Jupiters". 
When the star evolved on the giant branch it would become big enough to capture its 
planet/brown dwarf. The planet/brown dwarf spiraled into the envelope of the giant to 
initiate a CE phase. The frictional drag arising from its motion through the CE would lead 
to loss of its orbital angular momentum and deposit of orbital energy in the envelope. If 
there was enough orbital energy, the spiral-in proce ss would expel the envelope of the giant, 
leaving a WD remnant ( Nelemans fc Tauris 1998 1^. If the initial separation between the 



In addition, Bear fc Sokel ( 2010 ) suggest that the binary systems may reach stable synchronized orbits 
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star and the substellar object (s) is less than tens of Solar radii, the final outcome would be 
an under- massive WD with or without the surrounding planet /brown dwarf, depending on 
whether it evaporated, filled its own Roche-lobe, or survived. The discovery of a 0.053Mf7^ 



brown dwarf in a short (0.08 d) period orbit around a 0.39Mq WD 0137—349 (IMaxted et al. 
20061 ) presents strong observational evidence for this interaction. 



During the planet/brown dwarf- involved CE phase, the companion's envelope expanded 
rapidly and filled its Roche-lobe, leading to mass transfer onto the NS. Since this phase was 
very short (< 10^ yr), the mass transfer rate would be much higher than the Eddington 
limit rate for the NS, so that the NS accreted very small mass ~ 1O~^M0 (unless the CE 
phase lasted much longer time). If this is the case, we have to require that some MSPs were 
born this way, rather recy cled during the LMXB ev olution. In the literature, this idea has 
already been discussed by lMiller fc HamiltonI (120011 ) . who showed that, the existence of the 
innermost, moon-sized planet in the PSR 1257-1-12 system suggests that the pulsar was born 
with approximately its current spin frequency and magnetic field. A schematic view of the 
formation of BMSPs with planet/brown dwarf-involved CE evolution is shown in Fig. 8. 

Not only having an impact on the formation of MSPs, the substellar objects, if really 
exist in low-mass binaries, may also play a role in the CE evolution , espec ially infiuence 
the estimate of the CE efficiency parameter acE- Recently iDavis et al.l (120121 ) reconstructed 
the CE phase for the current sample of post-CE binaries (PCEBs) with observationally 
determined component masses and orbital periods. Searching for correlations between acE 
and the binary parameters, they found that, when the internal energy of the progenitor 
primary envelope is taken into account, acE dec reases with increasing mass Mp of the primary 
(i.e., the progenitor of the WD) (see however iDe Marco et al.ll201ll ). and c tcE ^, 1 for 1 5, 
Mp/M c7) < 2, which seems to be in contrast with uqe ~ 0.25 derived by iRicker fc Taam 
(120121 ) from numerical simulations. If there are planets/brown dwarfs around these low- 
mass primaries, they can contribute extra orbital energy to help expel the primary's envelope 
during the RGB/AGB phase, making the CE efficiency parameter within the canonical range 
< acE < 1. 



before the onset of the CE phase. Such stable synchronized orbits allow the RGB star to lose mass prior 
to the onset of the CE phase. Even after the secondary enters the giant envelope, the rotational velocity is 
high enough to cause an enhanced mass-loss rate. 
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4.2. The Rebirth Periods of BMSPs 



Accretion onto the NS in I/LMXBs changes both the mass and spin of the NS. The spin 
evolution of an accreting NS depends on the i nteraction between the magnetosphere and 



the accretion disk, as described as follows (e.g.. iBhattacharya fc van den Heuvellll99ll ): the 



accretion disk is truncated at the magnetospheric radius i?n 
radius of the disk 



which is close to the Alfven 



\2/7 



where i?s and Rs are the surface magnetic field and the radius of the NS, respectively. 
Material is channeled to the NS at the inner radius, producing a spin-up torque. If the 
NS is spinning very rapidly, the disk may be truncated outside the corotation radius R^o = 
(GMiP'^ /Ati'^Y^^ , and the NS exper iences a centrifugal barrier that can inhibit accretion 
(i.e., the so-called "propeller" effect; llUarionov &: Sunyaevlll975[ ). So the NS will eventually 



reach the equilibrium spin period such that -Rm = -Rco, or 



Peq ^ 2ABfRf''M;''/\MjM^^a)-'/'ms 



5/7/ 



(10) 



where Bg = BjlO^ G, and Re = Rs/10^ cm. This period can be regarded as the beginning 
or rebirth period of recycled pulsars when accretion terminates. 

Th e rebirth period of a MSP can be derived only when its actual spin-down time is 
known. iBhalerao &: Kulkarnil ( 120111 ) recently reported the optical discovery of the companion 
to the ~ 2Mq MSP J1614-2230. The opti cal colors show that the 0.5M ^ companion is a 2.2 
Gyr old CO WD. From the age of the WD. lBhalerao fc Kulkarnil ( 1201 ll ) calculated the period 
of the pulsar at birth, and found that pulsar should be born with a spin close to its current 
value, implying that the final accretion rate was < 10~^MEdd- This value is two orders of 



magnitude s n aaller than the estimate from theoretical calculations by iLin et al.l (1201 if ) and 
Tauris et al.l (120111 ). These authors suggested that the system began as an IMXB consisting 
of a NS and a ~ main-sequence secondary, which evolved to be a CO WD with He 

envelope. The NS gained the most mass during the final LMXB phase lasting ~ 5 — 10 Myr 
at near-Eddington rates. 

We present a systematic view on the relation between P^q/ Bg^"^ and Porb derived from 
theory and observations in Fig. 9. The solid lines outline the theoretically expected dis- 
tribution of Peq/Bg^'^ (a '^^'^ here Rq is taken to be 1) from binary evolution 
calculations, and the symbols repre sent the observations of recycled pulsars (data are taken 
from the ATNF Pulsar Catalogue; [Manchester et al.ll2005l ). Comparison between observa- 



tions and theory shows that, except in a few cases, the observational PspinZ-Bg^'' is generally 
larger than Peq/Pg^^ by ~ 1 — 2 orders of magnitude, challenging the simple recycling the- 
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oryQ. Thus, some other mechanisms must work to produce efficient spin down torque(s) or 
reduce the spin- up torque due to mass accretion, if the stellar evolution models are correct. 
The proposed explanations can be summarized as follows. 

(1) There is no definite spin equilibrium si nce GR may remove the ang ular momentum 
from the NS. This idea was ffist suggested by iPapaloizou fc Pringld (119781 ) to account for 
the cutoff in the spi n distribution of NSs in LMXBs. The ma in emission rn echanisms involve 
crustal mountains (IBildstenl 119981 ) . magnetic deformations (jCutler II2002I ). and unstable r- 
modes (lAnderssonI 119981 ) in the NSs. All these processes can produce a substantial mass 
quadrupole moment and thus a spin-down torque due to GR. The problem with this inter- 
pretation is that, recent observational results on some millisecond X-ray pulsars have shown 



Haskell & Patruno 


2011; 


Patruno et al. 


2012 



references therein). So the following models focus on modifications of the equilibrium period 
(Eq. [10]) that the NS will finally reach. 

(2) If the NS magnetic field lines can thread the a ccretion disk, an extra magnetic 
torque will be exerted on the NS ( Ghosh fc LambI 1979a ). and the equ i librium period wi ll 
be ^ 1 - 3 time s that in Eq. (10) JGhosh fc LambI Il979bl : IWaneJll995l : \u fc Wanejll996h . 
Andersson et al.l (120051 ) further suggested that the inner disk region may be geometrically 
thick and sub-Keplerian, and dominated by radiation pressure, if the mass transfer rate is 
above a few percent Eddington accretion rate. The coupling between the disk and magnetic 
field can reduce the amount of angular momentum deposited onto the NS from accretion 
by a factor A (the ratio of orbital angular velocity to Keplerian rotational speed in the 
disk, A < 1), compared with the c ase of thin disks {A = 1), thus a decreased spin-up 
torque results. iPatruno et al.l (120121 ) show that the existence of spin equilibrium as set by 
the disk— magnetosphere interaction is able to explain the observations of millisecond X- 
ray pulsars, if the spin-down torque coming from the interaction between the disk and field 
outside corotation is sufficiently large. Howe ver, the radial extent of the coupling between the 



star and the disk is highly controversial (e.g. iMatt fc Pudritzll2005l: lGhoshll2007l'). Stellar and 



disk winds may also take away the angular monientum of the NS (iGoodson et al.lll997l . Il999 
Ustyugova et al.l l2006l : iRomanova et al.l l2009l : IZhang fc Lill2010l ). In addition, the abrupt 



torque reversals obse rved in disk-fed X-ray pulsars Her X-1, Gen X-3, and 4U 1626—67 
(IBildsten et al.l 119971 ) suggest that current models of disk— magnetosphere interaction may 
have severe limitations. 



^Here we assume that the current spin periods of MSPs are not far fr om their initial ones . This holds if 
their characteristic ages are longer than the cooling ages of the WDs (e.g. lLorimer et al.lll995l) or the Hubble 
time. 
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(3) If Eq. (10) does apply, a decaying mass transfer rate at the end of the mass transfer 
process will result in relatively large rebirth period of BMS Ps, provided that the evo- 



lution time is longer than the NS spin evolution timescale (jjeffreyl Il986l : iRuderman et al. 



19891 ). More recently iTaurid (120121 ) suggested that this equilibrium will even be broken when 
IM2I decreases substantially, so that the magnetospheric radius i?m becomes larger than the 
corotation radius i?co, causing the NS to enter the propeller phase. A centrifugal barrier 
arises to expel matter entering the magnetosphere with a braking torque 



Nc^M2RinK{Rn 



(where QK{Rm) is the Keplerian angular velocity at R^) to act to slow down th e NS. An 
alternatively appropriate form for the accretion torque is (e.g. iMenou et al.lll999l ) 



N ~ -M^RlinKiRra) - a 



(12) 



where Qg is the angular velocity of the NS. With Eqs. (11) and (12) we recalcul ate the 
spin e volution of the NS when the donor star gradually decouples from its RL as in iTauris 
f l2012h . In Fig. 10 the blue and red lines denote the results calculated with Eqs. (11) and 
(12), respectiveljlf]. It is seen that the decaying mass transfer rate indeed causes an increase 
in the spin period, which is larger than the equilibrium period corresponding to the average 
mass transfer rate before decoupling. However, The calculated discrepancy seems not to be 
large enough to explain the deviation of the rebirth periods of BMSPs from that given by 
Eq. (10), especially if we adopt the torque form of Eq. (12). 

(4) A decaying magnetic field may also change the rebirth period of the BMSPs, if 
the field decay timescale is shorter than the spin evolution timescale, so that the final spin 
period deviates from the equilibrium period. During the decay of the NS magnetic field 
induced by accretion, the spin-up timescale increases as R^ decreases because the accreting 
matter carries less specific angular momentum, while the field decay timescale decreases as 
Rra decreases (fi eld decay stops when Rm = -Rs)i resulting in a s ignificant departure from 
the spin-up line ( iBurderi et al.lll996t iKonar fc Bhattacharyal 119991 ). 



(5) The rebirth periods of BMSPs may be related to the long-term spin equilibrium 
with varying mass accretion rate. A significant fractio n of LMXBs are likely to be transient 
systems subject to disk instability (jCoriat et al.l 120121 ). For example, an NS LMXB in a 10 
hr orbit would be transient if the average mass transfer rate is lower than ~ 10^^ gs^^. In 
addition, in narrow LMXBs, X-ray irradiation of the donor star could destabilize the mass 



■^The discontinuity at the point where the two hnes meets is due to the fact that we use Eq. (12) to 
calculate the spin evolution even before the propeller phase in both cases. 
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transfer, and lead to mass transfer cycles (iHameury et al.lll993[ ) - mass transfer is spasmodic 
with phases of high mass transfer driven by the thermal expansion of the convective envelope 
of the irradiated donor alternating with phases with low or no mass transfer, during which the 
donor readjusts towards thermal equilibrium of the un-irradiated star. The final spin period 
is thus determined by the spin- up during outbursts/high mass transfer phase balanced by 
the spin-down during quiescence/ low mass tra nsfer phase, which may be consider ably longer 
than that attained via stable mass accretion (ILi et al.lll998t iPatruno et al.ll2012l ). 



Currently our knowledge about the evolution of LMXBs and accretion disks is not 
sufficient to tell which one is the dominant factor in determining the rebirth periods of 
BMSPs. Perhaps they are influenced by a combination of (at least some of) the above- 
mentioned mechanisms. 



4.3. Origin of Intermediate-Mass BMSPs with Porb ^ 1 d 



It is well known that a large fraction of IMXBs will produce binary pulsars with a 
CO WD companion. As shown in Fig. 3, the orbital periods of these intermediate-mass 
binary pulsars (IMBPs) are generally >~ 3 — 6 d. This leaves the puzzle that the binary 
pulsars with Porb 5, 1 d and a CO WD companion are hard to explain reasonably in any 



binary population (jDeloyell2008l ). In Table 1 we list the parameters of IMBPs with P^ 



orb 



< 



1 d. Except PSR B0655+64, other pulsars share the characteristics of rapid spin and low 
magnetic field. 

The small orbital periods combined with massiv e WDs suggest that this type o f systems 
may have evolved through a CE and spiral- in phase (Ivan den Heuvel fc Taanj|l984l ). Similar 
to that of the Double NSs, their direct progenitors could be binaries consisting of a He star 
and a NS star in a close and circular orbit, as a result of spiral- in of a wide binary in late 
Case B or Case C mass transfer. The core of the companion was a He star which could 
already have gone through quite some He burning. After the spiral-in one then had a He 
star with already CO in its core. During He-shell burning, the envelopes of the He stars, if 
their masses are < 3.5 M©, will expand and slowly transfer mass to the NS, which will have 
caused the spin- up and recycling of the NS. Th is interpretation seem to be responsible for 
PSR B0655-I-64 (Ivan den Heuvel fc Taamlll984l ). which has the longest spin period (195.6 
ms) and highest magnetic field (P ~ 10^° G) among the pulsars in Table 1, implying that 
it has been partially recycled with little mass accreted. For other pulsars, the short periods 
and low magnetic fields (P ~ 1 0^ G) requir e that at least some 10~^ material has been 
accreted by the NSs. Recently IChen et al.l (l201ll ) presented the detailed binary evolution 
calculation for a binary consisting of a NS (of mass 1.3 Mq) and a low-mass He star (of mass 
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1.0 Mq) with an initial orbital period of 0.5 d. They showed that the mass transfer seems to 
be able to spin up the NS's spin to milliseconds, producing BMSPs like PSR J1802— 2124, 
which has a short orbital period Porb ^ 1 d. 

Another possible way to form compact IMBPs invokes evolution of IMXBs with anoma- 
lous MB. Traditionally MB is thought to work in low-mass, main-sequence stars. However, 
There are many intermediate-mass stars which have anomalously high magnetic fields (> 
100 G), i.e., Ap/Bp stars. In IMXBs, the irradiation-driven stellar wind from the surface 
of the Ap/Bp donor star by the N Ss may couple with t he high magnetic fields, resulting in 
angular momentum loss at a rate (jJustham et al.ll2006l ). 



MB 



,il)M2M-_ 



Dl5 

■'gmI' 



a/4 



(13) 



where B2 is the magnetic field strength at the surface of the donor star, and if) is a, parameter 
with typical value ^/c^ ^ 10~^, combi n ing th e wind-driving energy conversion efficiency and 



irradiation geometry. lJustham et al.l ( l2006l ) explored the effect of anomalous MB in the 
formation of compact BH LMXBs, in which the orbital energy of the companion star, if of 
low-mass initially, would be un able to expel the envelope of the B H progenitor. For the black 
hole LMXB XTE J1118-I-480, lOonzalez Hernandez et all (120121 ) measured a decay rate of 
the orbital period Porb = —1.83 ± 0.66 msyr~^. This value is much larger than predicted by 
most conventional MB and mass loss models for the binary parameters of XTE J1118-I-480, 
unless the companion star has a surface magnetic field P2 > 10 — 20 kG to enhance MB. 

In Fig. 11 we show the evolution of a binary initially consisting of a 1.8 Mq NS and a 3.5 
Mq donor with Porb — 2.5 d. We assume that the donor possesses a surface magnetic field 
B2 = 500 G. RLOF starts at the age ~ 2.4 x 10^ yr with a mass transfer rate ~ 1O~^M0 yr~^. 
Because the donor is initially more massive than the NS, the orbital period decreases rapidly. 
Most of the transferred mass is assumed to be ejected out of the binary in the vicinity of the 
NS, taking away both mass and orbital angular momentum. The orbit starts to widen when 
M2 ~ 2 Mq. When the donor mass decreases to be ~ O.8M0, the mass transfer rate reduces 
to be < 10"'^ M0 yr~^. In the subsequent ~ 8 x 10^ yr, MB begins to dominate the evolution. 
The orbital period keeps decreasing, the NS accretes most mass in this phase and is spun 
up to milliseconds. The final outcome is a tight binary system consisting of a ~ 2.1 Mq NS 
and a ~ 0.5Mq WD with Po^b ~ 0.4 d. 

The anomalous MB model will not work if the NS is initially of low-mass (~ IMq), 
since the allowed parameter space for forming binary pulsars will be very small (see Fig. 1). 
Additionally we note that models involving accretion-in duced collapse of an ONeMg WD 
cannot account for the Porb distribution of these BMSPs (ISutantyo fc Lill2000l ). 
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5. Summciry 

This work is motivated by the fact the evolution of I/LMXBs with canonical NSs seems 
to meet difficulties in explaining some of the observational characteristics of BMSPs, and 
the measurements of the NS masses indicate a wide distribution ~ 1 — 1.8 Mq. We have 
preformed numerical calculations of the evolution of I/LMXBs consisting of a 1.0 or 1.8Mq 
NS and a 1.0 — 6.0Mq donor star, to investigate its dependence on the initial NS mass, and 
the properties of the descendent BMSPs. The main results can be summarized as follows. 

1. The allowed parameter space in the initial Porb — -^2 diagram for forming recycled 
pulsars increases with increasing NS mass. This may help explain the formation of BMSPs 
with Porb ^ 60 d and their distribution in the Porb — diagram. Alternatively, some 
of these wide binary pulsars may be formed through mass transfer driven by planet/brown 
dwarf-involved CE evolution. 

2. The equilibrium spin periods Peq of accreting NSs in LMXBs derived from the 
standard magnetosphere- accretion disk interaction model are in general shorter than the 
observed spin periods of BMSPs by more than one order of magnitude. This implies that 
either the simple equilibrium spin model does not apply for the spin evolution in accreting 
NSs, or there are other mechanisms/processes to spin down the NSs when forming BMSPs. 

3. Some of the compact IMBPs might have evolved from IMXBs in which the companion 
star were strongly magnetized Ap/Bp stars with enhanced MB. 

4. Our calculations doubt the suggestion that the orbital period gap (~ 20 — 60 d) 
of BMSPs is related to the occurrence of the bump-related detached phase in the LMXB 
evolution, since the accretion disk would become unstable at earlier time. 

We are grateful to Ed van den Heuvel for helpful discussions and the referee for con- 
structive comments. This work was supported by the Natural Science Foundation of China 
under grant number 11133001, the National Basic Research Program of China (973 Program 
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Fig. 1. — The thick hnes confine the allowed parameter space in the initial orbital period 
vs. donor mass plane for I/LMXBs to successfully form binary pulsars with He and CO WD 
companions. In the left and right panels the initial mass of the NS is taken to be 1.0 
and 1.8 Mq, respectively. 
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Fig. 2. — Evolution of the NS mass Mi, the donor mass M2, the obital period Porb, and 
the mass transfer rate M2 for a IMXB consisting of a 1.0 Mq (upper tow panels) or 1.8 Mq 
(lower two panels) NS. 
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Fig. 3. — The final orbital period function of the WD mass Mwd for binary pulsars 

with an initial 1.0 Mq (left) or 1.8 Mq (right) NS. The number next to each curve denotes 
the initial mass of the donor star (the progenitor of the WD). The gray and black curves 
represent binary pulsars with He and CO WD companions, respectively. 
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Fig. 4. — Evolution of the NS mass Mi, the donor mass M2, the obital period Porb, and the 
mass transfer rate for a LMXB with a 1.0 Mq NS and a 1.0 Mq initial donor (upper two 
panels) or with a 1.8 Mq NS and a 1.5 Mq initial donor (lower two panels). The thick lines 
represent the critical mass transfer rate for unstable accretion disks. 




Fig. 5. — This figure shows the effect of outfiow from the L2 point on the formation of MSP 
J1614— 2230. The upper two and lower two panels correspond to 5 = and 0.04, respectively. 
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Fig. 6. — The orbital period as a function of the WD mass for MSPs with possible He WD 
companions. In the left and right panels the thick lines are theoretically expected relations 
for NSs with an initial mass of 1.0 Mq and 1.8 Mq, and the dots correspond to binary pulsars 
with an assumed mass of 1.2 Mq and 2.0 Mq, respectively. The error bars of the WD masses 
cover the 90% probability mass range for randomly oriented orbits having i = 90° toi = 26°. 
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Fig. 7. — The mass transferred rate (red dashed curves) and accreted masses (blue dotted 
curves) of the NS as a function of the final orbital period for binaries with a 1.0 Mq (left 
panel) or 1.8 Mq (right panel) NS. The black dot-dashed line show the critical accretion 
rate MEdd- In the left panel, the red dashed (blue dotted) curves, from bottom (top) to top 
(bottom), correspond to 1.0, 2.0, and 3.0 Mq donor stars, respectively. In the right panel, 
plotted are those in the cases of 2.0 and 3.0 Mq donor stars. 
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Fig. 8. — This schematic plot shows the possible formation process of wide BMSPs with a 
He WD companion. We assume that in a wide binary (with orbital period at least longer 
than tens of days) the secondary star initially possesses a (or multiple) substellar object (s), 
which spirals into the stellar envelope when it expands and develops a He core. The CE 
evolution may also cause the secondary to fill its Roche-lobe, and transfer mass rapidly to 
the NS. After the envelope is expelled, a wide binary with a NS and a He WD is produced. 
If the NS can accrete sufficient mass from its companion, it will appear as a MSP. 
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Fig. 9. — The black lines show the relationship between PeqZ-Bg^^ and Porb for all systems 
from our calculations. Blue diamonds, green stars, and magenta triangles represent binary 
pulsars with Pspin < 10 ms, 10 ms < Pspin < 100 ms, and Pspin > 100 ms, respectively. The 
red dashed lines from top to bottom correspond to mass accretion rate of 0.01, 0.1, and 
l.OMEdd for a IAMq NS, respectively. 
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Fig. 1 0. — Spin evolution of an accreting NS at the end of mass transfer phase as in iTauris 
(I2OI2I ). The black line denotes the evolution of the mass transfer rate, and the blue and 
red lines denote the spin evolution governed by Eqs. (11) and (12), respectively. The two 
vertical dashed lines confine the propeller phase. 
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Fig. 11. — The evolution of the NS mass, the donor mass, the orbital period, and the mass 
transfer rate for a 1.8 Mq NS IMXB. With anomalous MB, the final outcome is a IMBP 
with orbital less than 1 d. The thick line corresponds to the critical mass transfer rate for 
unstable accretion disks. 
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Table 1: Parameters of Iiitennediate-Mass Binary Pulsars with P^,^■l, < 1 d' 



Pulsar 


Pspin (ms) 


Porb (d) 




B{G) 


B0655+64 


195.6 


1.0287 


0.80 


1.17 X 10^° 


J1435~6100 


9.3 


1.3549 


1.08 


4.84 X 10^ 


J1757-5322 


8.9 


0.4533 


0.67 


4.89 X 10^ 


J1802-2124 


12.6 


0.6989 


0.98 


9.69 X 10^ 


J1952+2630 


20.7 


0.3919 


1.13 





"Data are taken from the ATNF pulsar catalogue. 
''Median mass for the inclination angle i = 60°. 



